Correlation Properties of anisotropic XY Model with a Sudden Quench 
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Starting from a general Hamiltonian which may undergo a quantum phase transition (QPT) with 
the change of a controllable parameter, we obtain a general conclusion that in a sudden quench 
system, when the final Hamiltonian is fixed, the behavior of the time-averaged expectation of any 
observable has close relationship with the gapless excitation of the initial Hamiltonian. To clarify 
our conclusion, we investigate the two-spin correlation of a XY chain in a transverse field under 
a sudden quench at zero temperature. The critical property of the derivative of quench two-spin 
correlation and the long-range correlation of the quench system are analyzed. 
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I. INTRODUCTION 

Quantum phase transition (QPT), which is driven by 
quantum fluctuation and occurs at zero temperature, is 
a very important research area recently [T]. Besides the 
static properties of QPT, the non-equilibrium dynamics 
induced by the quench of the parameter in the Hamilto- 
nian through a critical point has always been an attrac- 
tive topic in condensed matter physics [2TfT5] . 

Generally speaking, a dynamical evolution can be in- 
duced by either a sudden quench or a slow quench. Tak- 
ing the spin chain in the magnetic field as an example, 
the typical paradigm of sudden quench is as follows. Ini- 
tially (t < 0) the magnetic field h(t < 0) = a. Then at 
time t = the magnetic field is changed to b suddenly, 
namely h(t > 0) = 6 and the system begins to evolve 
from the initial state. After a long enough time evolution, 
the time-averaged expectation value of an observable A 
reaches a steady value, which we define as quench A. 
Quench A is (A(a, b)) = lim^^ | f*(A(a, b, r))dr [13], 
with (A(a, b, r)) the expectation value of observable A at 
time t = t with initial magnetic field a and final mag- 
netic field b. One can fix b(a) to investigate the relation 
between quench A and a (b). 

The magnetization M and its derivative with respect 
to magnetic field d^M of transverse Ising model are an- 
alyzed in Rcf. [T3]. It is found that dhM has similar 
behavior for static system and sudden quench system. 
Here interesting questions arise: why does this similarity 
happen and whether other quantities still exhibit such 
similarity? 

In this work, we start from a general Hamiltonian H (A) 
containing a controllable parameter A to search for the 
answers. We obtain a general conclusion that in a sud- 
den quench system, when the final Hamiltonian is fixed, 
the behavior of the time-averaged expectation of any ob- 
servable has close relationship with the gapless excitation 
of the initial Hamiltonian. There is a similarity between 
the critical phenomena of sudden quench and static sys- 
tem. Then we numerically investigate the sudden quench 
properties of two-spin correlation of an anisotropic XY 
chain in a transverse field [19] . The results are consistent 
with our conclusion. We also research the long-range 



correlation of the system in both directions parallel and 
perpendicular to the magnetic field. There is no quali- 
tative difference between static and sudden quench case 
in the direction parallel to the magnetic field, but a re- 
markable difference exists in the direction perpendicular 
to the magnetic field, which is consistent with the result 
in Ref.[3]. At last, we briefly discuss the case in which 
initial Hamiltonian is fixed. 



II. GENERAL MODEL 

Here we consider a general Hamiltonian H(X) = Hq + 
XHi which contains a controllable parameter A. We sup- 
pose the ground state of H (A) is non-degenerate and with 
the varying of A the system can undergo a QPT hap- 
pening at A = A c , where the energy gap of H(X) van- 
ishes in the thermodynamics limit (TL). This is the case 
for many condensed matter systems, especially for spin 
chains. For example H] — ^\ of and A is the magnetic 
field for transverse Ising model. 

Now we try to derive the singular behavior of the time- 
averaged expectation value of any observable A when 
QPT happens, both for static case and sudden quench 
case. At time t < 0, A = a and the eigenvectors and 
eigenvalues of H{a) are denoted as \ip n (a)) and E n (a) 
respectively. At zero temperature the system is in the 
ground state of this initial Hamiltonian H(a), namely 
\ty(t < 0)) = \ipo(a)). Therefore for static case, which 
means A = a for t > 0, the the time-averaged expec- 
tation value is just (A(a)) = (^o(a)\A\ip (a)). By the 
perturbation theory, we know 
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straightforwardly, where $t(z) means the real part of a 
complex number z. One can expect that d a (A(a)) will 
behave singularly at a = A c because the energy gap van- 
ishes in TL. 

If at time t = 0, A is suddenly changed to b and A = b 
is kept for t > 0, the system will evolve under sudden 
quench. The eigenvectors and eigenvalues of the final 
Hamiltonian H (6) are denoted as \<p n (b)) and u n respec- 
tively. The time-evolving wave function of the system 



at t > is |*(t > 0)) =En(Mb)\Ma))e- iUnt \Mb)), 
where we omit the relative phase. The time-averaged 
expectation value is 

1 f* 
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So we have 
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m^O n 



(4> n (b)\A\Mb)) 
E (a) - E m (a) 



xK[( J ff / ) Om (^ m (a)|0„(&))( < /»„(&)|V'o(a))]- 



(2) 



When b ^ A c , the numerator in the above Eq.Q is 
not zero. Therefore d a {A{a 1 b)) will show singularity at 
a = A c because of the vanishing of the energy gap in 
TL. However if b = A c = a, we have ip = (f>, so in 
Eq.Q (ip m (a)\(l) n (b)){(f) n (b)\ipo(a)) = 5 m , n 8 n ,o- Consid- 
ering m ^ in the sum, we can see both the numerator 
and denominator in Eq.Q are zero at a = b = A c , mean- 
ing that d a (A(a, A c )) maybe has no singular behavior at 
a = A c . 

Because Eqs.Q and ^ share the same denominator, 
we conclude that there must be similarity between the 
critical singular behaviors of d a (A(a)} and d a (A(a, b)}. 
The singularity and scaling behaviors come from the van- 
ishing of the energy gap of the initial Hamiltonian in TL. 
If there is a critical phenomenon of the initial Hamil- 
tonian with a vanishing energy gap which can be indi- 
cated by the singular behavior of d a (A(a)) at a = A c , 
d a (A(a,b)} will also have a similar singular behavior at 
a = A G if b 7^ A c . 



III. XY SPIN CHAIN MODEL 



There are two kinds of phase transition in anisotropic 
XY spin chain model. One belongs to the universal- 
ity class of the Ising phase transition and the other is 
anisotropic transition |10j . In this paper we focus on the 
Ising phase transition. The equilibrium Ising phase criti- 
cal points are \h\ = 1. InJ = l>0 case, the system has 
anti- ferromagnetic interaction with the staggered magne- 
tization in the x direction (a x ) s t ag = X^o^ - l)*(°f ) 
as the order parameter. For \h\ < 1, the system is in 
an ordered phase with \(a x } s ta S \ — 1, while for \h\ > 1, 
the system is in a disordered phase with (a x ) s ta,g = 0. 
Now we turn to the problem of dynamics under sudden 
quench. Using successive Jordan- Wigner, Fourier, and 
Bogoliubov transformations, the exact solutions of the 
evolved state of XY model under a sudden quench can 
be obtained |19) , by which the quench dynamical proper- 
ties of many observables such as magnetization and two- 
spin correlation can be studied. In our work we focus 
on the quench two-spin correlation, which is defined as 
G£(a,6) = lim^oo \ / (cr£(a, b, T)er£(a, 6, r))d,T, where 
a = x, y, z, n is the distance between two spins and a (b) 
is the initial (final) magnetic field. 



In this paper, we use the anisotropic XY model in a 
transverse field under a sudden quench at zero temper- 
ature to check our conclusion in the above section. The 
Hamiltonian of the spin- 1/2 anisotropic XY spin chain is 



JV-l 



H = E { 2 [v- + tW+i + Q- - > - ^'M}> 



where J is the coupling constant, erj, crj, crj are the Pauli 
operators at the jth lattice site and J(l ± j)/2 is the 
measure of the interaction strength in x(y) component 
between two nearest-neighbor spins. 7 is the anisotropy 
parameter which can change from (XX model) to 1 
(Ising model). h(t) is the time-dependent magnetic field. 
In our paper, we set J = 1 > for convenience and the 
periodic boundary condition cr^ = <7q , a = x, y, z is used. 
The phase diagram of this model is easy to be obtained. 



IV. QUENCH TWO-SPIN CORRELATION 

As the beginning we want to clarify our conclusion in 
Sec. II by considering d a Gi' (a, b), namely the nearest- 
neighbor two-spin correlation. Due to the difficulty in the 
analytical calculation of them, we choose to numerically 
solve them. We let 7 = 1 for convenience in this part and 
one should note that all conclusions below are also valid 
for 7 7^ 1. 

In FigfTJa), d a Gf(a,b 7^ 1) diverges at a = 1 for both 
the static and dynamical case. Moreover, we calculate 
that d a G\(a,b) — const x \a — with /1 s» 0.25, a 
universal exponent irrelevant with b (FigJIJb)). While 
as predicted in Sec. II, the sudden quench case with 
b = 1 is very special. We can see from FigJTJa) that 
d a G\(a, 1) is not divergent at a — 1. Similar to G\ (a, b), 



FIG. 1. 7 = 1. (a) d a Gl{a,b) for static case and dynamical 
cases b = 0,0.5,1,1.5. (b) ln(<9 a Gf (a, b)) as a function of 
In | a — 1 1 . One can see the slopes of the lines for different b 
are all the same, namely a universal constant -0.25. 




FIG. 2. 7=1. (a) d a Gi(a,b) for static case and dynamical 
cases b = 0, 0.5, 1, 1.5. (b) d a Gf (a, b) as a function of In \a— 1\. 
One can see the slopes of the lines for different b are all the 
same, namely a universal constant 




FIG. 3. 7 = 0.6. G z n (a,b) for (a) The static case b = a, (b) 
b = 0, (c) b = 1 and (d) b = 1.5. 



(a) Static 
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d a Gi(a,b 1) diverges at a = 1 for static and sud- 
den quench case. After numerical analysis, we can find 
d a Gf(a,b) = — In |a— 1|+/(&), with f(b) a constant which 
we do not concern with (Figj2]jb)). We also can see from 
Fig{2ja) that d a G x (a, 1) is not divergent at a = 1. 

Through the comparison between static and sudden 
quench cases of XY model, we find that in general sudden 
quench case (b ^ 1), the quench system and static system 
share similar critical scaling form of d a G z ^ x \a, b). This 
is consistent with our conclusion in Sec. II and reflects the 
memory of the quench system to the initial Hamiltonian. 

Now we shift our focus from the nearest-neighbor two- 
spin correlation to the correlation between two spins with 
arbitrary distance n. We show G^(a,b) as a function of 
a and n for both static and sudden quench case in Fig(3j 
where 7 = 0.6 is used. From this figure, we know as n 
increases, G^(a, b) reaches a constant quickly after a little 
increase of n. There is no qualitative difference between 
static and general dynamical case. 

For G^(a,b), there will be a big difference. We let 
7 = 1 for convenience here. In a static system (a = 6), 
from FigQa) we can see when a > 1, G x (a) tends to 
zero very quickly and is nonzero only if n is very small, 
demonstrating no long-range correlation in x direction. 
While in another phase a < 1, the ground state is degen- 
erate and possesses long-range correlation in x direction. 
However, in a sudden quench system, the case will be 



G"(a,b) 




(d) b=2. 



FIG. 4. 7=1. G^(a,b) in (a) static system and dynamical 
systems with (b) b = 0.5, (c) b = 1 and (d) b = 2. 



obviously different. If the final magnetic field b < 1, af- 
ter a long time evolution the long-range correlation in x 
direction only exists at the point a = b < 1 (that is just 
the static case), meaning that when b < 1 a change of 
magnetic field at t — will eliminate the long-range cor- 
relation in x direction. In FigQb), we plot the G^{a, b) 
up to n — 150 with b — 0.5. Through numerical fit of the 
two typical data G x n (0.4, 0.5) and G x n (0.45, 0.5) (Fig[5j), 
we find G^(a, b) ~ n~ Cl e~ C2 ™ , decaying to zero with n 
even when a and b is slightly different. So we think the 
G£(a, b) is destroyed after an arbitrarily small quench. 
While b > 1, no matter what a is, even for a < 1, there 
is no long-range correlation in x direction. This means 
that there is a possibility that a system possessing a long- 
range correlation in x direction at the beginning will lose 
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FIG. 5. Numerical fitting of G*{a,b) at b = 0.5 up to 
n — 150. Through numerical fit of the two typical data 
G£(0.4, 0.5) and G£(0.45, 0.5), we find G x n {a, b) ~ n- cl e - C2n , 

With G£ (0.4, 0.5) = n -0.0000404236 e (-0.074916n-0.0031414) and 

G a: (0 45 5) = n - - 0000460715 g(-°- 07259538n - ' 00080686 ) \Y e 
extrapolate those two curves to n = 400. 



it after a long enough time evolution (FigQd)). So, only 
when static case a = b < 1 , the long-range correlation in 
x direction is preserved. In other cases we can have when 
n — > oo, (o-q (a, b)a^(a, b)) — > 0, leading to (a x (a,b)) = 
straightforwardly. A physical intuition tells us that be- 
cause the magnetic field is in the z direction, the long- 



range correlation in the z direction can survive after the 
quench but the long-range correlation in the x direction 
is destroyed by the quench. 

We note that the long range correlation is also re- 
searched by Sengupta et al. in Ref.[3]. They give some 
analytic results of the two point correlation function per- 
pendicular to the magnetic field direction in quantum 
Ising model. However, they only consider two limit cases, 
namely the initial magnetic field is fixed at a = and 
a = oo. In both cases, G„(a,b) tends to zero when 
ri — > oo but in different ways which depend on the value 
of the final magnetic field b. In Fig|4j one can see that 
(0,0.5) tends to zero much slower than G^(0, 2) and 
there is a clear spatial oscillation of G*(2, 0.5). These 
phenomena are consistent with Ref.[5]. Therefore, our 
numerical results confirm their analytic results about 
long range correlation and give an extension to the gen- 
eral situation beyond the two limit cases. 



V. A FURTHER DISCUSSION: FIXED INITIAL 
MAGNETIC FIELD 

Up to now we only consider the case in which the fi- 
nal magnetic field b is fixed and all quench quantities 
are regarded as functions of initial magnetic field a. For 
completeness, we discuss the case in which a is fixed and 
all quench quantities are regarded as functions of b. As 
one can expect, the quench quantities will have different 
behaviors. If we make a similar calculation with that in 
Sec. II using perturbation theory, we can obtain 



d b (A(a,b))=2j2J2 
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It is difficult to study the behavior of Eq. ^ near b = \ c 
for a general a. However, the case in which a = X c is 
easy to analyze. When a = A c = b, <j) = ip. There- 
fore in the first term of the right hand side of Eq.(|3]), 
{ipo{a)\(l>n{b)}((l>m(b)\i>o(a)) = S„ h0 5 nfi , meaning that 
this term is because m =/= n in the sum. Similarly, 
we have (ipQ(a)\(j) n (b))(4> n (b)\tpo(a)) — S n .o in the second 
term of Eq. (131) . So we have 



d b (A(\ c ,b))\ b = Xc =2Y, 



m^0 



9K[(ffj)om(<MA c )M0 o (A c ))] 
w (A c ) - oj m (X c ) 



This equation is the same with Eq.|I]) and means that 
db(A(X c ,b)) will diverge at b — A c because the vanish of 
the energy gap. 

In the following we use 9(,Gf(a, b) as an example to 
analyze (see Figj6|. As we expect above, dbGf(a,b) is 



divergent at b = 1 when a = 1. While a ^ f, it is 
discontinuous at b = 1. How to explain the discontinuity 
is still an open question to us. 



VI. SUMMARY 

In summary, starting from a general Hamiltonian 
which may undergo a QPT with the change of a con- 
trollable parameter, we obtain a general conclusion that 
in a sudden quench system, when the final Hamiltonian 
is fixed, the behavior of the time-averaged expectation 
of any observable has close relationship with the gapless 
excitation of the initial Hamiltonian. This general con- 
clusion, clarified by our example of XY model, to a large 
extent explains the similarity between the critical phc- 
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FIG. 6. 7 = 0.6. dbGl(a,b) with fixed initial magnetic field 
a = 0,0.5, 1,2. 



nomena of sudden quench and static system. The long- 
range correlation is also studied and we find that sudden 
quench can destroy long-range correlation in the x direc- 
tion. As we know, the one-spin reduced density matrix 
can always be written as p = (I + J2 a =x y z ( cr a) (T a)/^- If 
considering symmetry breaking effects in XY chain [20] . 
we cannot suppose (a x ) — 0. Therefore for static case 
it will be a complex task to solve (a x ). But for sudden 
quench, (er x (a, b)} = for a b. So our result can help to 
simplify the calculation of one (two)-site reduced density 
matrix. This will be useful in studying the entanglement 
of one site with others for XY model. In addition, the 
general conclusion can be useful for other quantities and 
other models. 

This work was supported by NSF of China un- 
der Grants No. 10821403 and No. 10974234, pro- 
grams of Chinese Academy of Sciences, 973 grant No. 
2010CB922904 and National Program for Basic Research 
of MOST. 



[1] S. Sachdev: Quantum Phase Transition, Cambridge Uni- 
versity Press, Cambridge 1999. 

[2] M. Greiner et al, Nature(London) 419, 51 (2002); T. 
W. B. Kibble, J. Phys. A 9, 1387 (1976); W. H. Zurek, 
Nature(London) 317, 505 (1985). 

[3] K. Sengupta, S. Powell and S. Sachdev, Phys. Rev. A 69, 
053616 (2004). 

[4] W. H. Zurek, U. Dorner and P. Zoller, Phys. Rev. Lett. 

95, 105701 (2005). 
[5] J. Dziarmaga, Phys. Rev. Lett. 95, 245701 (2005). 
[6] A. Polkovnikov, Phys. Rev. B 72, 161201(R) (2005). 
[7] A. Sen, U. Sen and M. Lewenstein, Phys. Rev. A 

72,052319 (2005). 
[8] B. Damski and W. H. Zurek, Phys. Rev. A 73, 063405 

(2006). 

[9] F. M. Cucchietti, B. Damski, J. Dziarmaga and W. H. 

Zurek, Phys. Rev. A 75, 023603 (2007). 
[10] V. Mukherjee, U. Divakaran, A. Dutta and D. Sen, Phys. 

Rev. B 76, 174303 (2007). 
[11] K. Sengupta, D. Sen and S. Mondal, Phys. Rev. Lett. 

100, 077204 (2008). 
[12] A. Polkovnikov and V. Gritsev, Nat. Phys. 4, 477 (2008). 
[13] M. Rigol, V. Dunjko and M. Olshanii, Nature (London) 

452, 854 (2008). 
[14] Y. Li, M. Huo and Z. Song, Phys. Rev. B 80, 054404 

(2009). 

[15] M. Rigol et al, Phys. Rev. Lett. 98, 050405 (2007); M. 

Rigol et al, Phys. Rev. A 74, 053616 (2006). 
[16] D. Patane, A. Silva, F. Sols and L. Amico, Phys. Rev. 

Lett. 102, 245701 (2009). 
[17] S. Hofierberth, I. Lesanovsky, B. Fischer, T. Schumm, J. 

Schmiedmayer, Nature 449, 324 (2007). 
[18] G. Roux, Phys. Rev. A 79, 021608(R) (2009). 
[19] E. Barouch, B. M. McCoy and M. Dresden, Phys. Rev. 

A 2, 1075 (1970); E. Barouch and B. M. McCoy, Phys. 

Rev. A 3, 786 (1971); E. Barouch and B. M. McCoy, 

Phys. Rev. A 3, 2137 (1971). 



[20] T. R. de Oliveira, G. Rigolin, M. C. de Oliveira and E. 
Miranda, Phys. Rev. A 77, 032325 (2008). 



